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The damping of collective modes and the relaxation of temperature anisotropics in a trapped Bose 
gas is determined at temperatures above the Bose-Einstein condensation temperature in the coUi- 
sionless regime. We demonstrate for both cases how the effects of collisions may be treated pertur- 
batively and calculate relaxation rates based on a variational principle. Our results are compared 
with experiment and with previous theoretical calculations. 
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I. INTRODUCTION 

In recent years a number of experiments have been car- 
ried out to study relaxation processes in trapped Bose 
gases. Such experiments may be divided into two broad 
categories, those on collective oscillations and those on 
temperature anisotropics. Relaxation rates of collective 
modes have been measured in rubidium vapors by Jin et 
al. §,|l, and in sodium ones by the MIT group |g]. The 
relaxation of temperature anisotropics has been investi- 
gated by Monroe et al. ^ and Arndt et al. g in cesium, 
Newbury et al. M] and Myatt Iq] in rubidium, and Davis 
et al. |9| in sodium. 

In the experiments on the relaxation of collective 
modes ||l|-^ , the atomic cloud is excited at the frequency 
of the normal mode of interest by means of a weak ex- 
ternal perturbation. The amplitude of the oscillation is 
then measured in the absence of the perturbation, and 
from that the damping rate of the mode is extracted. In 
the experiments on temperature relaxation, the atomic 
cloud is either prepared in a nonequilibrium state with 
different effective temperatures for motion along and per- 
pendicular to the symmetry axis of the trap l5l-S] , or it 
is driven out of equilibrium by displacing the trap center 

Most experiments on relaxation have been carried out 
in the coUisionless regime, where the collision rate is small 
compared with the frequencies of particle motion in the 
trap and also, in the case of collective oscillations, the 
oscillation frequency. In addition the mean free path is 
large compared to the size of the cloud. In this paper we 
focus on the coUisionless limit, at temperatures above the 
condensation temperature, Tc^ but our results also apply 
to the condensed phase at temperatures close to the tran- 
sition temperature, where the thermal excitations consti- 
tute the dominant component of the gas. A brief account 
of some of our results has been given previously [|l^ . 

Relaxation of temperature anisotropics in trapped 
clouds of alkali-atom vapors has been studied theoret- 
ically both analytically [^,0 and numerically |p|,p|,pl| . 
Newbury et al. and Myatt (s) used the Boltz- 



mann equation with an ansatz for the form of the 
non-equilibrium distribution function to obtain an an- 
alytic expression for the temperature relaxation rate and 
solved the problem numerically for large temperature 
anisotropics. Monroe et al. ||] performed a Monte-Carlo 
simulation of temperature relaxation for atoms with an 
energy-independent scattering cross section, and Arndt 
et al. made a similar calculation for resonant scattering 
g. Finally Wu and Foot solved the Boltzmann equa- 
tion numerically and extracted the rate of temperature 
relaxation pi[ . 

The present paper is based on the semi-classical Boltz- 
mann equation, from which we derive an eigenvalue equa- 
tion for the relaxation rate. Treating the collisions per- 
turbatively, we obtain general expressions for the damp- 
ing rate of the collective modes and the relaxation rate for 
temperature anisotropics. We include effects of quantum 
degeneracy and explore the role of energy dependence of 
the scattering cross section. 

The paper is organized as follows: In Sec. II we develop 
the general formalism. We then examine in Sec. Ill the 
damping of modes with frequencies equal to twice the 
characteristic frequencies of the trap and calculate the 
damping rate for a general, energy-dependent cross sec- 
tion. The case of an energy-independent scattering cross 
section is studied in detail, and analytical results are ob- 
tained in the limit of Boltzmann statistics for a simple 
trial function. We perform calculations with improved 
trial functions, and find that the corrections are small, 
of the order of 1%. Our results are then compared with 
experiment. 

In Sec. IV we examine the relaxation of temperature 
anisotropics and show that this problem is essentially the 
same as that of the damping of collective modes. Again 
we perform calculations with a simple trial function and 
with a more general one. The improved trial function 
yields a reduction of the relaxation rate by about 7% for 
the two cases which we examine. We then compare the 
results with data from experiment and from numerical 
simulations. Section V is a brief conclusion. 



II. GENERAL FORMALISM 

When the energy of a particle is much greater than 
the spacing between energy levels the motion may be de- 
scribed semiclassically. In a trap, energy spacings are of 
order the oscillator quantum, hcoi, where tUi is the fre- 
quency associated with motion in the i direction. Even 
when T is as low as the transition temperature, the ra- 
tio kBTc/huj ~ [Af/C(3)]"'^/'^ is very much greater that 
unity, where N is the total number of atoms in the trap, 
C{a) is the Riemann zeta function, and iD — {cuxUjyUJzY'^, 
and therefore the level spacing is small provided the 
anisotropy of the trap is not extreme. Also we may 
neglect the influence of the mean fleld due to interac- 
tions between atoms. This gives rise to a potential of 
order n(r)Uo, where n(r) is the number density and 
Uq = Anh a/m is the effective two-particle interaction 
with a being the scattering length for atom-atom colli- 
sions, and m being the atomic mass. Even when T = Tc, 
the energy nUo is typically only about 1% of ksT, and 
the molecular fleld may therefore be neglected in the 
equation of motion of the particles. 

Since the motion can be described semiclassically the 
state of the system may be specified in terms of a particle 
distribution function /(p, r, t), where p is the momentum 
vector, r is the position vector, and t is the time. This 
satisfies the Boltzmann equation. 



|[ + p . v/ - vy • Vp/ = 

at m 



I, 



(1) 



where V is the external potential, and / is the collision 
integral, which is a functional of /(r, p, t). We write the 
distribution function a,s f = f^ + 6f, where /" is the 
equilibrium distribution function and Sf is the deviation 
from equilibrium, and linearizing Eq. dlf) we find 

^f + ^-V5f-VV-VpSf = -m], (2) 

where $ is defined by Sf = /"(I -I- /°)<1>. 

We now consider collisions. The energy levels of a 
single particle moving in a harmonic trap are discrete, 
and for a spherically symmetric trap there is strong shell 
structure. For the anisotropic traps used in experiments, 
the shell structure is less pronounced, and the density of 
states will be made more uniform by interatomic inter- 
actions. The ratio nUo/hCu is of the order of unity for 
typical experimental conditions. Finally collisions will 
smear the density of states. As a consequence, we expect 
that under current experimental conditions the collision 
integral will be well approximated by the expression for 
a bulk system. The linearized collision integral is thus 

/[$] = / dn / da |v - vi|($ + $1 - $' - $'i) 

x/°/i°(i+r')(i+/i"'). (3) 

Here the momenta of the incoming particles in a colli- 
sion are denoted by p and pi, and those of the outgoing 



ones by p' and p[, and dr = dp/{2Trh)^. In the above 
expression we have introduced the differential cross sec- 
tion da, which in general depends on the relative velocity 
u = |v — Vi| of the two incoming particles, as well as on 
the angle between the relative velocities of the colliding 
particles before and after the collision. Since the atoms 
obey Bose-Einstein statistics, their distribution function 
in equilibrium, /'^(p,r), is 



r(p,r) = 



1 



p^/2m+V-fi]/kBT _ 2' 



(4) 



where fj, is the chemical potential. We shall consider har- 
monic traps where the potential has the form 



TA /'22i 22, 2 2\ 



(5) 



In the absence of collisions the motions of particles in 
the three coordinate directions decouple, and particles 
execute simple harmonic motion for each of the coordi- 
nates. Let us consider motion in the z direction. In the 
absence of collisions (/[<&] = 0) the general solution of 
the Boltzmann equation having time dependence 6""*"*=* 
has the form 



51 = [Pz + imcOzzy-hiKj), 



(6) 



where h is an arbitrary function, and Kj are constants 
of the motion in the absence of collisions. So, neglecting 
collisions, for Sf{r,p,t) = h{Kj){pz + imujzz)^e~^^'^'*, 
one has 



d 

— + niujz ) (5/ = 0. 



(7) 



We now consider the effect of collisions. Provided the 
collision rate is small compared with the oscillator fre- 
quency, collisions will not mix modes of the system in 
the absence of collisions with different values of n. Con- 
sequently we expect the solution to be of the form 



(5/(p,r,t) == h{K.j){pz+imujzzTe'' 



,t-Tt 



(8) 



This problem therefore amounts to doing degenerate per- 
turbation theory as in quantum mechanics, because we 
consider mixing of modes which have the same frequency 
in the absence of the perturbation, in this case collisions. 
Collisions will generally couple particle motions in differ- 
ent directions, but we begin by assuming that the colli- 
sion rate is small compared with \lOz — i^x\ and \ujz — Wy|, 
in which case this coupling may be neglected. It is im- 
portant in traps with axial symmetry for modes in which 
particles move perpendicular to the symmetry axis, and 
we shall consider its effects later. 

Substituting Eq. (^ into the Boltzmann equation one 
finds 



d_ 
dt 



df = -TSf = -m. 



(9) 



Equivalently, Eq. (^) may be written as 



(10) 



The problem thus reduces to solving the eigenvalue equa- 
tion (|lO|) and it can therefore be treated by conventional 
perturbation theory, with the collision integral being the 
perturbation. We are interested in the mode with the 
smallest decay rate, since this is what is measured exper- 
imentally, and therefore, rather than solving Eq. (HG) di- 
rectly, we use a variational method. A related variational 
method is known to give highly accurate results in the 
calculation of transport coefhcients |E2| . The most long- 
lived mode corresponds to the smallest value of Tj, which 
we denote by Fq. Using Eq. { \[(\ ) for any trial function ^, 
we multiply by ^* and integrate over the coordinates and 
momenta, and find that 



rn< 



(**/[*]) 



(l*p/°(l + /"))' 



(11) 



which is our basic variational expression for the damp- 
ing rate. Here the brackets denote integration over both 
coordinate and momentum space. The numerator in the 
above equation expresses essentially the entropy genera- 
tion rate, whereas the denominator gives the energy in 
the mode ^. 

It is worth pointing out the connection between 
Eq. (O) and the usual variational expression for trans- 
port coefficients that applies in the static limit. A trans- 
port coefficient such as the shear viscosity is determined 
by an inhomogeneous Boltzmann equation of the form 



f{l + f)X^I[E] 



(12) 



and the transport coefficient is proportional to the quan- 
tity (S*X/0(1 + /0))p, where 

{E*Xf{l + ,f))p= /dr/0(l + /0)S*(p)X(p). (13) 



The above integration is over momentum space alone, in 
contrast to the previous case [Eq. (pTJ)], which involves 
integration over both coordinate and momentum space. 
Then for any trial function VP one finds |12| that because 
the eigenvalues of the collision integral are never negative, 



{E*xfii+n) 



< 



{^*imi 



(vl/*X/"(l+/0))2- 



(14) 



Comparing Eqs. (|ll| ) and (|j), we see that both are 
given by the same sort of variational expression and in- 
volve the collision integral. The differences are in the 
normalization condition which appears in the denomina- 
tors of the two expressions, and in the spatial integrations 
inEq.(|nD. 



III. DAMPING OF MODES 



A. The damping rate 



The inequality (p_l|) may be used to calculate the damp- 
ing rate of any mode in which particle motion in the z 
direction is perturbed, but let us focus on a trap with uni- 



axial symmetry, ojx 



and on modes with frequency 



2llJz, which are among those of experimental interest. In 
the same spirit as one adopts in calculating transport co- 
efficients, we begin by making the simplest choice, Eq. (g) 
with n = 2, and with h put equal to unity. Thus as our 
initial trial function $ we use 



^ — {Pz + imoJzzY 



(15) 



In Sec. Ill C we examine more general trial functions. 



The function <i> contains terms of the form 



ZPz 



and 



pI. It is clear from Eq. (||) that the collision integral 
gives zero when operating on the first two of these terms 
because collisions conserve particle number and momen- 
tum; it also gives zero when operating on p^ because of 
the conservation of the kinetic energy. Therefore we may 
write 



{{p. 



'nipz 



']) 



{pli[pl])^{{pl~p'/3)i[pi~py3]). 



(16) 



The second of these forms follows as a consequence of 
kinetic energy conservation, which allows us to substract 
any multiple of p^ from p^ in /[p^]. The choice made 
above introduces the function p^ — p^/3, which is or- 
thogonal to the collision invariants, since pi — p^/3 — 
(2/3)p^P2(cos 6), where Pi{x) is the Legendre polynomial 
of order I, and 9 is the polar angle of the momentum vec- 
tor. Substitution of Eq. ( p^ into Eq. ( |ll|) leads to the 
result 



rn< 



{{pl-py3)I[pl-py3]) 
'[pl + (mc.,z)2]2/0(i + /O)) ■ 



(17) 



This may be expressed in terms of the viscous relax- 
ation time r^^var calculated within the variational approx- 
imation. To see this, we start with the expression which 
connects T^,var with the viscosity 77, 



1 



77=-(|A|2/"(l + /"))pr,,var, 



(18) 



where the trial function has been put equal to A, i.e., the 
driving term in the Boltzmann equation. The variational 
expression for the viscosity is 



V ^ < 77vai = "1 



(X*/[A])p 



{\x\^P{i + P))l 



(19) 



and therefore we find with use of Eqs. (nSl) and (nS) that 
Eq. (n^) can be written as 



"-6 /dr((p2_p2/3)2jo(i + /o))p • ^ "i 
The factor of 1/6 comes from the ratio of the integrals 

reflects the fact that collisions damp directly only the 
part of the distribution function varying as pi — p^/3, 
while all parts of the distribution contribute to the total 
energy associated with the oscillation. In the classical 
limit Eq. (ISO) takes the simple form 



r°^^ 



dr 



n(r) 



drn{r) 



(21) 



We return to Eq. (MTtj now, starting with the numerator 
on the right side of the inequality. This quantity includes 
integration over real space as well as momentum space, 
and we write it as 



{{pl-py3)i[pl-py3]) 



drlc[z{r)], 



(22) 



where z(r) = z(0)e~^/'^^"^ is the fugacity, with z(0) = 
gH/kBT j^gJQg j^g value at the center of the cloud. /c[z(r)] 
is given by 



/,[z(r)]= dr dTi da|v-vi| 

x$*(a> + $1 - $' - ci>;)/"/l'(l + /°')(1 + fi')- (23) 

In the above expression unprimed quantities refer to in- 
coming particles and primed ones to outgoing particles. 
Using the symmetry of the collision integral under inter- 
change of incoming and outgoing particle momenta, we 
may write Eq. (p3) as 



lMr)]^^JdTjdnJdaW-Vi\ 

X 1$ + $1 - $' _ $;|2/0/0(l + /"')(1 + ff). (24) 

Because of the Galilean invariance of the collision pro- 
cess it is convenient to use the center-of-mass momen- 
tum P and the relative momentum coordinates p^ and 
p^ instead of the momenta of the incoming and outgoing 
particles p,Pi, and p',p'i. 



P = p + 

Pr = 7:(P 



Pi =P 

-Pl), 



Pi; 

p; 



1 



2V- --' -'■ 2 
We introduce A[<I>] by the definition 



(P' - P'l) 



(25) 
(26) 



and for $ 



Pz 



p^/3 one finds 



Ab2_p2/3] = 2(p2,_p;/) 

1 . 



= 2 
4 



pI.z 



^Pr 



1 ,2 
oPr 



-pi[P2{cOs9r)-P2icOS0',)] 



(27) 



where 9r and 9'^. are the polar angles of the vectors Pr 
and p'j. with respect to the z axis. We write the element 
of cross section da in Eq. (^3|) as {da/dn,)dfl, where dO is 
the element of solid angle for the vector p^. Since we are 
considering scattering of indistinguishable bosons, the fi- 
nal states corresponding to pj, and — pj, are identical, 
and we perform the summation over final states by sum- 
ming over all Q and then dividing by 2 to avoid double 
counting. Using the addition theorem for Legendre poly- 
nomials and their orthogonality property, we find from 
Eq. (|27|) that the numerator of Eq. (p^) is 



((p^_pV3)/b^-pV3]> 



16 

A5m 



dr dr dri 



x^Jdn[l-P2{cosx)]^plff,a + f)i^ + fi'). (28) 

where x is the angle between p^ and pj,. The product 
of the Bose functions in Eq. ( P8|) depends on the angles 
between P and the vectors p^ and p'. , which implies that 
it is impossible to factorize the collision integral, except 
for the case of Boltzmann statistics. In order to perform 
the angular integrations we express sin x i^^ terms of 9^ 
and 9[,, using the addition formula 



cos X — cos t)r cos ( 



sin Or sin 9'j. cos{(j)r ~ 4>'^), (29) 

where 0^ and 0J, are the azimuthal angles of the vectors 
p,. and p'j.. Equation ([2^) then gives 



2ir 



di<Pr - K) Sin' X - ^(1 + l/ + y" - Sy^y'^), (30) 



with y = cos 9r and y' — cos 9'^. . 

The differential cross section da/dfl depends in general 
on y and y' and in principle wc may consider any form 
of da/d^. At the low energies of interest it is, however, 
a good approximation to assume that the scattering is 
purely s wave |§-|3jljl • When the cross section is energy- 
independent, which is a very good approxiation for ru- 
bidium and sodium atoms, da/dfl — 4a^. The factor of 
4 reflects the fact that for identical bosons the scatter- 
ing amplitude is twice that for distinguishable particles. 
For cesium atoms there is a bound state close to zero en- 
ergy and the scattering is resonant. If the resonant state 
were at zero energy, the s wave phase shift So would be 
TT (modulo 27r) and therefore the cross section at low en- 
ergy would be given by the unitarity limit for s waves, 
da/dO. — 4//c^. More generally, when the resonant state 
is close to, but not exactly at zero energy, one may take 
the leading term in the effective range expansion p5[ and 
write tan(5o = ka. Since da/dVt — (4/fc^)sin (5o, one 
then finds 



da 

dh 



4a" 



1 + {kascY 



(31) 



where k — Pr/h is the wavevector corresponding to Pr- 
The differential cross section thus depends on the en- 
ergy of the colliding particles. For s wave scattering the 



integral of the P2 term in Eq. (^8|) vanishes, and the in- 
tegrand on the right side in this equation involves the 
total momentum-dependent cross section (t{u), which is 
defined as 



a{u) 



dn 



da 

dn' 



(32) 



where for the relative velocity of the incoming particles, 
u, u = 2pr/in. For a differential cross section of the form 
of Eq. (bTI) the total cross section is 



a{u) 



Sna^ 



1 + (fca)2 



(33) 



We turn now to the denominator of Eq. (17) and define 
D as 

D = JdrJdT [pl + {mLo,zf]'fil + f) 

^2ldrldT [pt+pl{mi,,zr]f{l + f), (34) 

where the latter result follows from the fact that /° is 
a function of p^ + {mujzz)'^, which is symmetric under 
interchange of Pz and itilUzZ. The integral in Eq. ( |34| ) is 
over a six-dimensional hyperspace with coordinates p and 
mcur, and the "angular" and radial integrals decouple. 
With the use of the equation fil+f) = -ksTdf/dE, 
with E being the particle energy, one finds 



D 



dr dTp^f{l + f) 
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^mkeT I dr / drp'^f 



where 



2i5/2;r 



D{v) 



{mkBTy^ f drD{r), 



(35) 



X dx 



z-i(r)e^ 
-5-35/2 [2; (r)], 



(36) 



with5„(z) 



^j^ z jl^ being the familiar Bose integral 



From Eqs. (l2§), (|3^), (|33|) and (|35|) we obtain 
4 



rn< 



dr dr dri 



15m?kBT , 
xpla{u)ff,{l + r){l + f>-) / I dr I drp'f . (37) 



Our final variational inequality (p^) can be written in the 
form 

1 f mkBT Y^' _ 



TM^)]dr / / drD{r) , (38) 



which is one of our main results. Using Eq. (^6|) we 
perform the spatial integration in the denominator of 
Eq. (|3|) and find 



rn< 



1 



{^] TT^ Ic[z{r)]dr. (39) 



307:7/2^2 V h } 54 [^(0)] 



In the above equations ctq = Svra^ is the scattering cross 
section at zero energy, and 



V = (8fcBT/7rm)i/2 



(40) 



is the mean thermal velocity. The quantity /c is propor- 
tional to /c and it is given by 

Ic = dxr / dxa dy dy' x'^x].a{xr,t) 

Jo Jo Jo Jq 

,,2 I „/2 o„,2„/2n 



X 



F{xr, xo,y, y') (1 + y' + y" - ?>y^y'^). (41) 



In the above expression we have introduced a{xr, t) which 
is given by the ratio a{u)/ao and it is a function of u = 
Pra/h. We have also defined Xr = Pr/imkBTY^^ and 
t = T/Tq, with fcsTo = ^ /ma?' , which implies that u — 
Xri^^"^. Also xq is given by a;o = P/{imkBTy^^. Finally 
we have introduced the product F = f'^f i{l + /'^ )(1 + 
fi ), which may be written in the useful form 



F = 



y^ p (^0+^r) 



(1 — ze ^i)(l — ze ^2)(1 — ze ^3)(1 



(42) 



where Xi = pi/{2mkBT)^''^ with pi being the momentum 
of particle i. The indices i — 1 and 2 refer to incoming 
particles, and z = 3 and 4 to outgoing ones. The variables 
Xi are explicitly given in terms of Xr,xo, y, and y' by 

A = 2(^0 + 2x0X^2/ -I- xl), xl = -{xl - 2xQXry + xl), 

^3 = 2(^0 + '2xoXry' + xl),xl = -{xl - 2xoXry' + xl). 

Taking the classical limit of Eq. (^8|) , we find that the 
integrand in the numerator is proportional to z^(r) and 
therefore to the square of the density, as expected. The 
integrand in the spatial integration in the denominator 
is linear in the density in this limit. In addition, in the 
classical limit the integrations over the center of mass 
and the relative momenta decouple, and Eq. ( pq ) may be 
written in the form 



rn< 



30%/2 



(Jqv 



I drnl^ir) drn^iir) 



X / dxr x^e '"^ a{xr,t) 



dxx^e-"" , (43) 



where 7ici(r) = nci(0)e^^/'^^-'" is the density in the clas- 
sical limit, with nci(O) being the density of the atoms at 
the center of the cloud, where it is a maximium. 



To express Eq. (^) in a more physical form it is conve- 
nient to define an average density weighted by the parti- 
cle number, or 



nl{r) dv / / nei(r) dv = nei(0)/23/2 , (44) 



where the last equality holds for a harmonic potential. 
The result Eq. ( |43| ) may be written in a form that makes 
it physically more clear, 



1 



To < -^(r{u)u I drnliir) 



rfrrici(r) 



< 



24^2 



nci{0)a{u)u, 



where 



a{u)u 



/p duu"^ exp{—mu'^ /AkBT)(T{i 
L duu^ exp{—mu'^/4:kBT) 



(45) 



(46) 



Using Eq. (^5|) we find that classically for energy- 
independent s-wave scattering, 



To < Tr«cl(0)o-oU = To.cl. 

15 



(47) 



We turn now to the effects of quantum degeneracy. For 
energy-independent scattering and for T ^ Tc one may 
derive an analytical expression for the damping, as shown 
in Appendix A, and the result is 



1 , , _ 
To < — »^cl(0)cro« 
15 



^C(3) (I 



(48) 



where nci(O) = NCJ^[m/2'KkBT]^/^. The coefficient 
3C(3)/16 is approximately equal to 0.23. The dashed 
line in Fig. 1 shows the result of the expansion given by 
Eq. (|4^), as a function of Tc/T. Under other conditions 
the integrals in Eq. ( p8[ ) must be calculated numerically. 
The relaxation rate can be written as 



To 
Fo.ci 



= G{TJT). 



(49) 



The solid line in Fig. 1 shows the function G{Tc/T) for s 
wave, energy-independent scattering. The value of Fq 
at Tc is « 1.66Fo,ci. The assumption of an energy- 
independent cross section is a rather good approxima- 
tion for Rb and Na atoms, but later we shall consider 
the case of an energy-dependent cross section, which is 
relevant for Cs. 
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FIG. 1. Solid line: The relaxation rate Fq given by Eq. ( |3q ) 
for energy-independent scattering, as a function of Tc/T, nor- 
malized to the classical value, Fo.ci- Dashed line: The result 
of the analytic expansion given by Eq. ( ^Iq) also normalized to 

Fo.ci. 



B. Traps w^ith degenerate modes 

Up to now we have studied oscillations along one of 
the axes of the cloud, taken to be the z axis, and we have 
assumed that the mode is non-degenerate: ojz 7^ w^, and 
^z 7^ ^y There are experimental situations, however, 
where the damping of oscillations perpendicular to the 
axis of the trap has been measured W,01. Since Ux = ^y, 
collisions couple the motion in the x and y directions. Be- 
cause the collision integral is invariant under rotations of 
the spatial coordinates about the z direction, the appro- 
priate solutions of the Boltzmann equation may be char- 
acterized by the "magnetic" quantum number, ttIz, spec- 
ifying the symmetry under rotations about the z axis. 
From the trial functions (px+imujxx)'^ and (py+imajyi/)'^ 
for motion in the x and y directions one can construct 
two sorts of trial function 

$(m2 = 0) = {px + imuixxY -f {py + imuyy)^, (50) 

and 

^(Im^l = 2) = {px + iniLUxx)'^ - {py + imujyyY . (51) 



For TUz = the only part of this trial function that con- 
tributes to the collision integral is p^. + Py, and since 
Px+Py = P^ ^ pI: the collision integral gives the same 
result as in the case studied in Sect. Ill A. On the other 
hand with this choice of trial function the denominator 
is equal to 2D, where D is given by Eq. (p^), because the 
energy in the mode is twice as large. This is easily seen 



by calculating {^{mz = 0)p and performing the angular 
integrations, which yield 

JdvJdTmmz^O)\'fil + f) = 

2 fdv fdr [pl + im^^xf]'fil + f) 



dr dTp^f{l + f), (52) 
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since the integral of the cross terms vanishes. The de- 
nominator in the present problem is thus twice as large 
as the one obtained in Eq. (pq), and one finds 



r(m, -o) = iro. 



This result is independent of the degree of quantum de- 
generacy. 

Turning to the trial function with jm^l = 2, Eq. (51) 



we see that it is a superposition of terms with niz — ±2. 
Since the ruz = +2 and iriz = —2 modes are degenerate, 
this mode has the same damping. The only part of the 
trial function that contributes to the collision integral is 



Px 



Py . Making use of the expressions 



Px 



' -l^{8n/15)'/V{Yi+Y,-') 



(54) 



and 



pI+pI- 2pV3 = (2/3)V2(8Vl5)VVr2^ (55) 

where YJ™ are the spherical harmonics, we find for the 
damping rate r(|mz| = 2) of the \mz\ — 2 mode that 



ro(|m,|=2) = 3ro(TO, =0) 



(56) 



since the denominator for <i>(|?Ti2| = 2) is the same as for 
^{niz = 0). The above ratios are also independent of the 
quantum degeneracy. 



C. Improved estimates 

One can improve the estimates of mode damping by us- 
ing more general trial functions. We start with Eq. (|^), 
which gives the most general form of Sf. In the calcula- 
tion presented in Sec. Ill A we assumed that the function 
h{Kj) = 1, or <i>o = {pz + imojzz)'^- More general trial 
functions can be used, and one systematic method for 
generating these is to expand h in terms of complete sets 
of functions of the variables Kj . For the problem under 
consideration the conserved quantities are the energies 
Ei ~ pf/2m + muJlrl 12 for each of the coordinate axes. 
For a trap with an axis of symmetry the angular momen- 
tum about this axis is conserved, but this will not be a 
relevant variable, since in the experiments there was no 



preferred direction of rotation about the axis. The num- 
ber of possible functions of Ei is vast, and we shall not 
perform a general search. Rather we shall look for im- 
proved functions that depend only on the total energy. 
Also for simplicity we consider the classical limit only. 
We expand h in terms of a set of orthogonal polynomi- 
als, and a convenient choice is the Soninc polynomials 
introduced in the theory of transport coefficients in clas- 
sical gases ll^ . Because of the form of the normalizing 
factors in the denominator in Eq. (O) the appropriate 
orthogonality condition for the polynomials is 



{h,h,E^f) = J, 



jjj 



(57) 



where the factor E'^ reflects the presence of the prefactor 
(53) p'2 _|_ (777,^^2;)^ in the trial function, and Stj is the Kro- 
necker delta function. The functions hi are thus propor- 



tional to the Sonine polynomials 5| |1§| . We remark that 
the normalization we have chosen, which is convenient for 
our purposes, differs from the one used for Sonine poly- 
nomials. The simplest example is /iq = 1, and the next 
one is hi — \/5(l — E/5kT) and thus 



$1 = V5{pz+imLUzz)^{l- E/bkeT). 



(58) 



As our improved trial function we take a linear combina- 
tion of ho and hi, which is equivalent to c -I- -E, where c 
is a constant. We write the improved trial function as 



$i, 



(1 



)^/^$o + a$i. 



(59) 



where a is a variational parameter. We proceed by cal- 
culating the damping rate Fo.imp corresponding to this 
trial function and minimizing it with respect to a. In 
terms of the quantities 



7,,, = (**/[*,]) 



jj/' 



(60) 



we obtain 



ro,imp = ]^ [(1 - "')7o,o + 2a(l - a^y/^ji^o 

+a^7i,i] ■ 



(61) 



For the calculation of the quantities 71^0 and 71.1, we no- 
tice that the only difficult part is the one that involves 
the combination {pz + imoj zz)^ E / ^ksT . One can eas- 
ily see that due to the conservation laws, only the parts 
{p^ 12711 + V)p^ + iujzzpzp^ give non-vanishing terms, in 
addition to the terms we considered earlier with <i>o only. 
We find that in terms of the coordinates Pr, pj-, and P, 
the combination of deviation functions (26) is given by 



A[$i„p] = 2 (p2 ^ - p; /) [(1 - a2)i/2 ^ ^^ 



a 



V5 V aV5 



1 



5 ksT 10 rnksT 
oi2^/l 1 ( P, 



P2 



5 ksT \2m 



lOJzz] {pr,z'P-Vr-p'rz'P-V'r). (62) 



which reproduces Eq. (^^ for a — 0. We calculate the 
matrix elements and find that 71,0/70,0 = — v5/20, and 
71,1/70,0 = 199/80. Identifying the ratio jo,o/D as the 
quantity To which we calculated earlier, we write Eq. ( |6l| ) 
as 



To, imp (a) = To 



7o,o 7o,o 



1 ^ n 2N1/2 . 119 2 
1 a 1 — ct I ' H a 

10 ^ ' 80 



(63) 



Minimizing the above expression with respect to a, we 
find that a « 0.073; the minimum is « 0.992 To and 
therefore the correction is of the order of 0.8%. Since 
a is very small, it is a good approximation to expand 
Eq. (p3) in powers of a and write To. imp in the form of a 
perturbative expansion as in quantum mechanics: 



- 0,imp 



7?,o 



70,0(71,1 - 7o,o) 



(64) 



The correction term 1 — 7^ 0/70,0(71,1 — 7o,o) is equal to 
118/119 « 0.992. Our final' result is therefore 



To,. 



118 
1785 



nci(0)o-ow- 



(65) 



This provides evidence for the result of the earlier calcu- 
lations being very close to the actual value. 



D. Comparison with experiment 

The MIT group has reported two measurements of the 
damping rate of oscillations along the axis of a cigar- 
shaped cloud. In the first one Q iV = 5 x 10^ atoms 
[ p9[ were observed to oscillate at a frequency of « 35 Hz 
at T = 2 Tc and the damping time was measured to be 
about 80 ms. The frequencies of the trapping potential, 
which was axially symmetric, were v^ ~ 19 Hz along the 
axis of the trap, and i>±^ w 250 Hz perpendicular to it. 
Using these numbers we find that classically the central 
density is w 2.9 x 10^^ cm"'^ and assu ming that a — 27 .h 
A ||lj] for Na atoms, wc find from Eqs. ( p9| ) and ( pSJ ) that 
the damping time in the coUisionless regime is ~ 46 ms. 
However, the system is neither in the coUisionless nor in 
the hydrodynamic regime. The oscillation frequency is 
« 35 Hz, which is intermediate between the one in the 
coUisionless limit vc — 2vz ^ 38 H z, and the one in the 
hydrodynamic limit vh ~ ^Jl2/b v^ w X.hhv^, « 29 Hz 
pO| . We can use the simple interpolation formula given in 
Ref. ||2l| to calculate the corrected damping time. Ac- 
cording to this formula, the damping time is approxi- 
mately (yc — vii)l(y — V}i) times the value we found 
earlier in the coUisionless limit, which gives a damping 
time of « 70 ms, in rather good agreement with the ex- 
perimental value of sa 80 ms. As shown in Refs. [E2[, 



a detailed study of the frequency and the damping of 
the modes in the regime intermediate between the hy- 
drodynamic and the coUisionless limits shows that for the 
modes studied experimentally the interpolation formula 
we introduced in Ref. pit] is a very good approximation. 

In the other experiment performed by the MIT group, 
the number of atoms was iV « 8 x 10^ at T « Tc and 
the frequencies of the trapping potential were v^ k, 230 
Hz, and v^ k, \1 Hz. The oscillation frequency was 
measured to be v k, 1.75 i^z and the damping time to 
be K, 50 ms. The peak density calculated classically is 
nci(O) « 6.3 X 10^'^ cm~'^ and our estimate Eqs. ( ^ ) and 
(pq) for the damping time for this experiment in the col- 
lisionless limit is w 18 ms at T = Tc- Using again the 
interpolation formula from Ref. |2^ we find that the cor- 
rected damping time is ~ 41 ms, which is also in reason- 
able agreement with the experimental value k, 50 ms. 

The JILA group has measured the damping rate of 
oscillations in the plane perpendicular to the axis of the 
trap H. The trap frequencies in this experiment are 
cja; = tijj^ = cjz/v8 = 27ri^_L, where v\^ « 129 Hz. For 
iV = 8 X lO'' at T = l.lTc we find that classically the 
central density is « 3.7 x 10"'^^ cm~^. Assuming that 
a = 53 A for rubidium atoms JE3] , our calculations of 
the relaxation times Eqs. (^9|), (|53D, ( p6| ) and (^) imply 
that the rriz = mode has a damping time of ~ 88 ms, 
and the \'mz\ = 2 mode has a damping time of « 29 
ms. The observed damping time was, however, w 50 
ms for the rriz — mode, and ~ 100 ms for the \raz\ = 2 
mode, indicating the need for further work to resolve this 
discrepancy. 



IV. RELAXATION OF THERMAL 
ANISOTROPIES 

A. Form of the trial function 

We turn now to temperature relaxation. As we men- 
tioned in the Introduction, in the experiments described 
in Refs. pHS] the energies associated with the motion 
of atoms in the various directions in the trap are un- 
equal and to a first approximation these anisotropics can 
be characterized by ascribing different temperatures to 
the motions in different directions. This problem is es- 
sentially the same as the one we presented earlier. The 
most general form of the change in the particle distri- 
bution function Sfx is in this case any function p of the 
energies Ei along the direction i, 



5fT^p[E.,,Ey,Ez). 



(66) 



It can be checked easily that this form of SJt satisfies 
the Boltzmann equation in the absence of collisions with 
a frequency equal to zero. If one changes the temperature 
associated with particle motion in the z direction by an 
amount 5T and that associated with motion in the trans- 
verse direction by an amount — (5T/2, conserving in this 



way the average particle energy and the density up to 
first order in ST, the change in the particle distribution 
function S/t is thus 



SfT cx /0(1 + f) 



Ez - -x{Ex + Ey) 



«^/°(l + /°) 



E 

Ez 

3 



(67) 



Since the results do not depend on the anisotropy of the 
trap, we assume that the external potential V is isotropic 
with a frequency m. As our initial trial function we 
therefore take a function proportional to E-, — E/i — 
(pI - pV3)/2m + mtu^z^ - r^/3)/2, 



<^T ^ P^ P2{9p) + imcor)^ P2{9r), 



(68) 



where 9p is the polar angle of the momentum vector, and 
6r is the polar angle of the position vector. Later we shall 
consider improved trial functions. 



B. Calculation of the thermal relaxation rate 

Let us now calculate the thermal relaxation rate Ft- 
For the trial function ^t given by Eq. (psh the only 
part of it that gives a non-zero contribution when acted 
on by the collision integral is the term proportional 
to p^cos^(6'p) = pI- Since in the previous problem 
of the damping of oscillation modes we assumed that 
$ = (p^ + imUzzY, the quantity ($^ J[$t]) for ther- 
mal relaxation is the same as that for the relaxation of 
oscillations, Eq. (|2q). 

Turning to the denominator Dt = (|$Tp/°(l + f)) 
which appears in the thermal relaxation rate, the cross 
terms that arise from squaring $t vanish because of the 
orthogonality of the Legendre polynomials in Eq. (68), 
and therefore 



Dn 



Ft < 



Ft < 



dv j dT [p^ + {mu:rY]f{l + f) 
dvjdrp^fil + f). (69) 

D/5. 

(70) 
(71) 
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Comparing the above equation with Eq. (p5f) , Dt 
Finally we find that 



i'^Tmr]) 



3Fn 



2 /dr((p2_p2/3)2j-0(i + /0))p • 



This result may be understood by noting that the de- 
cay rate for temperature anisotropics in a uniform gas 
is given by Ft < t~1^^ in the corresponding approxima- 
tion. In a trap there are two facts that must be taken into 



account. The first is that the system is inhomogeneous, 
and therefore the relaxation time must be averaged over 
space. This accounts for the ratio of integrals in Eq. ([Tl]). 
The second is that the total energy associated with a tem- 
perature anisotropy has contributions from the potential 
energy as well as from the kinetic energy. For an har- 
monic trap, the average potential and kinetic energies 
associated with motion along a given axis are equal and 
therefore the total energy for a given direction is twice the 
kinetic energy. Since collisions isotropize directly only 
the kinetic energy, the time to isotropize both the kinetic 
and potential energies is twice as large as that calculated 
neglecting the effect of the potential energy. This is the 
origin of the factor of 1/2 in Eq. ([7l|). 

Neglecting the effects of quantum degeneracy, Eq. (|71 
gives 



Ft < 



1 



8^/2 



nci{0)(T{u)u. 



(72) 



Equation ([72[) gives for an energy-independent cross sec- 
tion cr(w) — (To, 



Tt < -nci{0)aov. 
5 



(73) 



Considering the case of a more general form for the scat- 
tering cross section, Eq. ([72[) can be written in the form 



(74) 



rT<^nci{0)aovF{t). 
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As defined earlier, t — T/Tq, with fcsTo = h /ma? being 
the characteristic energy that enters the scattering cross 
section and 



F{t) = 



dx. 



7 - 



(75) 



For the specific choice a = {\ + x1t) ^ , the above expres- 
sion takes the form 



F{t) = 



l + xH 



(76) 



For T <^Tq the scattering cross section is essentially en- 
ergy independent, a{u) ~ (Tqi whereas in the opposite 
limit the cross section is resonant, (j{u) ~ %-k jk^ . For Cs 
atoms, a lower bound on the magnitude of their scatter- 
ing length has been determined to be 260 as Q, where 
as is the Bohr radius; therefore Tg < 15 //K. Figure 
2 shows the result of Eq. ( [74| ) for Cs atoms, for a peak 
density nci(O) = 10"'^'' cm~ . The x axis shows the tem- 
perature in units of Tq and the y axis shows Ft measured 
in units of 



F 



T,0 



1 

15 



nci(0)cro«O, 



(77) 



where Wq is the velocity wq — {ST:kBTQ/m) 
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FIG. 2. The thermal relaxation rate Ft as a function of 
T/To, as given by the classical version of Eq. (k4), for a cross 
section a{k) = 87ra^/[l + (fea)^]. Ft is measured in units of 
rT,o = nci(0)aowo/15. 




FIG. 3. The effective cross section (Toff given by Eq. (plb, 
as a function of T/Tq, for a cross section of the form 
o'(k) — Sira^/ll + (ka)^]. The cross section is measured in 
units of (Tcff.o given by Eq. (|82|). 



A simple interpolation formula for the ratio Tt/^t,o is 
given by 



3\/t 



Ft.o [1 + (3f) 



ivll/l/ ' 



(78) 



where ly ~ 0.9. For Cs atoms we find that Ft,o > 0-2 s^^, 
assuming that the magnitude of the scattering length is 
greater than, or equal to 260 a^ Q and for nci(O) = 10^° 
cm~'^, which is the typical density in the experiments 
reported in Ref. Q performed with Cs atoms. 

Equation (|7^ ) for Ft is useful for the case of Cs atoms; 
under the experimental conditions H it is a rather good 
approximation to assume that the scattering is resonant, 
with a — Sir/k'^, where k = mu/2h is the relative mo- 
mentum, since T > Tq. We then obtain the result in the 
classical regime 



^ 64 nei(O)^' 
It < -7 j^- 

15 TO^U 



(79) 



It is instructive to define an effective cross section Ccff 
according to 



fcff 



Tt 



(80) 



From Eq. (74) we see that 



Ceff 



2V2 

IE' 



aoF{t). 



(81) 



Figure 3 shows a^g given by Eq. ( p(\ ) as function of T /Tq. 
For energy-independent scattering, 



Ccff 



2\/2 



■O-Q — fcff.O, 



whereas for resonant scattering. 



Ccff 



128%/2 h^ 



15 (rnu)^ 



(82) 



(83) 



Finally using Eq. (bTI) we plot in Fig. 4 the effective cross 
section CTcff for Cs atoms as a function of the tempera- 
ture for four different values of the scattering length, a — 
100 as, 200 as, 300 as, and 400 as- For very small values 
of the temperature tXcff = (7cfi,Q — 0.142, 0.569, 1.278 and 
2.272 xlO~^^ cm^, respectively. 
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FIG. 4. The effective cross section acs for Cs atoms 
given by Eq. (pll), as a function of T, for a cross section 
of the form a{k) = 87ra'^/[l + (ka)'^]. The four lines cor- 
respond to the following choices for the scattering length, 
a = 100 as, 200 as, 300 as, and 400 as from bottom to top. 



C. Improved estimates 

Following the same method as in Sec. Ill C we cal- 
culate the corrections to the thermal relaxation rate Tt 
using better trial functions. We choose the improved trial 
function to be 



(1 - a^)^/^$T,o + a^TA, 



(84) 



where a is a variational parameter. ^t,o is given by 
Eq. ( |6^ ) and for $t,i we take 

*T.i = \/5 [{pI - pV3 + (mujfiz^ - rV3)] 

X (1 - E/5kBT) (85) 

by analogy with Eq. (jSg) . In this case we find that only 
the parts proportional to Ep^—p'^/6m give non- vanishing 
terms in addition to the terms we considered earlier with 
^T.o only, and thus 



A[$T,imp] = 3(p: 
aV5 V 



r,z 






y')'/'+aV5 



5 kgT 



p2 

4 



10 mkBT 

(p.p,,.p . p,. - p.p;,^p ■ p;.) 

For the case of energy-independent scattering, we find 
that the matrix elements 71,0/70,0 = — v5/20, and 



10 raksT 

a^/5 1 



71,1/70,0 = 277/240. Therefore for energy-independent 
scattering the thermal relaxation rate Frimp correspond- 
ing to $T,imp is given by 



rT,imp(a) = Ft 



v5 ., 2\l/2 , 37 2 

10 ^ ' 240 



(87) 



The minimum is at a ~ 0.47 and the value of the function 
is « 0.941 Ft, i.e., the correction is « 6%. Therefore, our 
final answer for FT,imp for the case of energy-independent 
scattering is 



T.imp 



0.189 nci(0)croi^. 



(88) 



In a study of the problem of temperature relaxation, My- 
att 1^ has used the Boltzmann equation to obtain an 
analytic expression for the temperature relaxation rate, 
and has solved the problem numerically for large temper- 
ature anisotropies. The coefficient 1/5 which appears in 
Eq. (^) is consistent with his results. In a numerical ap- 
proach of the same problem, Monroe et al. [^j have found 
that a Monte Carlo calculation gives for this coefficient 
a value of 0.185, which is remarkably close to our result, 
0.189. 

For the case of resonant scattering we find 71,0/70,0 = 
\/5/20, and 71,1/70,0 = 53/48 and thus the correction is 
« 7% in this case, or 



nci(O)^^ 



(89) 



As shown by Arndt et al. in Ref. |6[|, a Monte Carlo 
calculation for F^ for the case of resonant scattering gives 
a value for this coefficient of « 64/(2 x 10.7) w 64/21 w 3, 
instead of the coefficient 4 that appears in Eq. (p9) . 



D. Comparison with experiment 

The experiment performed by Arndt et al. [|6| with Cs 
atoms shows that the collisional cross section is inversely 
proportional to the temperature in the range 5 — 60 /iK. 
This indicates that To ^ 5 /xK. Since these temperatures 
are larger than Tq for Cs, according to our study the cross 
section should indeed vary as T~^ . Our result given by 
Eq. (H) fits the data of Ref. § rather well. 

In the experiment of Monroe et al. [B the collisional 
cross section of Cs atoms was measured to be w 1.5 x 
10""'^^ cm^, and it was independent of the temperature in 
the range between « 30 — 250 /iK. This value of the cross 
section is of the same order as the data in Ref. g for the 
highest temperatures. As Fig. 4 implies, our calculated 
cross sections for these temperatures are also of order 
10-12 cm2. 



In the experiment of Davis et al. |9|] , a quadrupole trap 
which gives a potential proportional to {x^ +y'^ +4:Z^y^'^ 
was used to confine sodium atoms. The atoms were ob- 
served to oscillate with an effective frequency of the or- 
der of kHz. A non thermal distribution was created by 
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temporarily displacing the trap center along the z axis. 
The time for relaxation to equilibrium was measured and 
found to be of the order of seconds, for a temperature of 
200 /iK and for densities in the center of the trap of order 
lO^*' cm-3. 

We cannot immediately apply our results to this exper- 
iment because the potential is not separable in cartesian 
coordinates. However Monte Carlo simulations show that 
motions are effectively separable (non ergodic) for peri- 
ods of at least 10 s Q. To estimate the relaxation time 
for temperature anisotropies, we observe that the rate 
at which the particle distribution is isotropized in mo- 
mentum space by collisions is given by the ratio of the 
integrals in Eq. (^Tj) . If the potential were separable in 
cartesian coorinates and were a homogeneous function of 
the coordinates of degree v^ the potential energy associ- 
ated with motion in any direction would be v jl times the 
kinetic energy, and therefore the relaxation rate allowing 
for the potential energy would be 



T,v 



< 



1 



Jdrr-XA{pl-pV3ff{^ + f)) 



l + iy/2 Jdr{{pl-py3rp{l + n)p 



(90) 



This result reduces to Eq. ( [7l|) for an harmonic poten- 
tial {v = 2), and to the expression Tt < Tnl^r for a 
uniform gas {v — 0). Because the mixing times due to 
the non-separability of the trapping potential are long, 
we would expect that for the quadrupole trap this result, 
with V = 1, should give a good estimate of the thermal re- 
laxation rate. For classical statistics, energy-independent 
scattering and v — I, Eq. (pO[) gives 



Tta< 



15V2 



nc\{0)aov. 



(91) 



The spatial average that occurs in the present case is 
/ n2j(r) dr/ / nci(r) dv = nci(0)/23, instead of ?ici(0)/23/2 
for harmonic potentials. 

The relaxation rate that was measured in Ref. can 
be written as 



rate of relaxation processes in a homogeneous gas is that 
the distribution function has a more complicated spatial 
variation. This is because the motion is essentially col- 
lisionless, and consequently the distribution function is 
a function of both momentum and position, whereas for 
a homogeneous gas in the hydrodynamic limit the rele- 
vant distribution function is a function only of momen- 
tum, apart from a rather simple spatial dependence. We 
formulated the problem as a variational one, and found 
that for the simplest trial functions the damping times 
for vibrational modes and for thermal relaxation differ 
only by numerical factors. With the simplest variational 
functions, the most long-lived mode in thermal relaxation 
corresponds to an anisotropy of the temperature, as Mon- 
roe et al. assumed in their calculation. This picture is not 
true in general, and an improved trial function for this 
mode led to a damping time 6% longer than the simplest 
variational estimate, and in agreement with what was 
found in Monte Carlo simulations of the coUisional dy- 
namics for a constant cross section. For a resonant cross 
section there appears to be a significant discrepancy be- 
tween our results and the Monte Carlo simulations of the 
coUisional dynamics of Ref. Q, and the reason for this 
is at present unclear. For damping of vibrational modes, 
the simplest trial function gives a very good description, 
and an improved one gave changes in the damping time of 
less than one per cent. If necessary, further improvement 
of our estimates is possible by employing more general 
trial functions. The existence of rather precise expres- 
sions for relaxation rates should prove a useful tool in 
analysis of experimental data. 

Our calculations were carried out under the assump- 
tion that the semiclassical approximation is valid, not 
only for the streaming terms in the Boltzmann equation, 
but also for the collision term. For smaller numbers of 
particles, and/or larger trap frequencies the effects of the 
discreteness of the spectrum of single-particle energies in 
the trap will become important, and the situation will be 
closer to that encountered in studies of collective motion 
in finite nuclei 1241. 



T,l 



■nci{0)crov, 



(92) 



if we assume that a — 27.5 A for sodium atoms p4[ . 
Comparing Eqs. ( pi| ) and (JO^) we see that there is a dif- 
ference of a factor w 4.3. We regard the agreement as 
satisfactory because of the uncertainty in our estimate 
due to the unknown role of the lack of separability of 
the potential, and because of the difhculty of measuring 
particle densities experimentally p9|. 



APPENDIX A: QUANTUM DEGENERACY 
CORRECTIONS 

Including effects of quantum degeneracy to leading or- 
der in the fugacity z(r), we find for energy- independent 
scattering, cr(w) — ctq, that the collision integral in 
Eq. (|3|) gives 



V. CONCLUDING REMARKS 

In this paper we have calculated the damping of low 
vibrational modes and of thermal anisotropies using a 
Boltzmann equation approach. The main difference be- 
tween the present problem and that of calculating the 



Qz{r)]^z'V^h+z^^^\. (Al) 

Expanding the denominator of Eq. ( |38| ) we find 

D(r) =z^(l + z^] . (A2) 
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Combining Eqs. (Al) and (A2) we thus find that Eq. 
takes the form 



ro< 



rnksT 
15\/2 \2^ 



3/2 



(JQ W X 



^^(r) l + .(r)A^^Ur 




(A3) 



Performing the spatial integrations in Eq. (A3), we find 
for a harmonic potential 



To < TT cTou — -— 2- 
15 V 27r?i^ 



3/2 



z(0) 



5z(0)' 
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(A4) 



where z(0) is the fugacity calculated at the center of the 
cloud, r = 0. We now calculate z(0) as a function of the 
total number of particles A^ — J2r, I '^^ f'^i^'P)-. which 
gives up to order 2:^(0) 



N 



keT 

Tiuj 



m 1 



.(0) 



(A5) 



Solving the above equation in terms of z{Q) up to N"^ and 
using the expression for the critical temperature ksT^ = 
to[iV/C(3)]i/3, we find 



z(0) 



C(3) 



C(3) 



Inserting Eq. (A6) into Eq. (A3) we obtain 



To < -j-rnci(0)crow 
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C(3) 



(A6) 



(A7) 



where nci(O) = Nuj^[m/2'KkBTf^'^ ■ 
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